Dynamical fidelity of a solid-state quantum computation 



(N 
O 

o 

C 
3 



(N 



> 
00 

in 

o 
o 

^ : 

-i— > ■ 
G ■ 
cd 

3 : 
cr 



13 



G.P. BermanW, F. Borgonovil 6 ^, G. CelardoM, F.M. Izrailev^, and DTKamenev^ 

' a l Theoretical Division and CNLS, 
Los Alamos National Laboratory, 
Los Alamos, New Mexico 87545 
^ Dipartimento di Matematica e Fisica, 
Universita Cattolica, via Musei 41, 25121 Brescia, Italy 
[c] I.N.F.M., Unitd di Brescia and I.N.F.N., 
Sezione di Pavia, Italy 
ld] I.N.F.M., Unitd di Milano, Italy 
' e ' Instituto de Fisica, Universidad Autonoma de Puebla, 
Apdo. Postal J-48, Puebla 72570, Mexico 

In this paper we analyze the dynamics in a spin-model of quantum computer. Main attention is 
paid to the dynamical fidelity (associated with dynamical errors) of an algorithm that allows to create 
an entangled state for remote qubits. We show that in the regime of selective resonant excitations of 
qubits there is no any danger of quantum chaos. Moreover, in this regime a modified perturbation 
theory gives an adequate description of the dynamics of the system. Our approach allows to explicitly 
describe all peculiarities of the evolution of the system under time-dependent pulses corresponding 
to a quantum protocol. Specifically, we analyze, both analytically and numerically, how the fidelity 
decreases in dependence on the model parameters. 
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I. INTRODUCTION 

Many suggestions for an experimental realization of 
quantum computers are related to two-level systems 
(qubits). One of the serious problems in this field is a 
destructive influence of different kinds of errors that may 
be dangerous for the stability of quantum computation 
protocols. In the first line, one should refer to finite tem- 
perature effects and interaction of qubits with an environ- 
ment Jj|. However, even in the case when these features 
can be neglected, errors can be generated by the dynam- 
ics itself. This "dynamical noise" can not be avoided 
since the interaction between qubits and with external 
fields are both necessary for the implementation of any 
quantum protocol. On the other hand, the inter-qubit 
interaction may cause the errors. Therefore, it is impor- 
tant to know to what extent the interaction effects may 
be dangerous for quantum computation. 

As is known from the theory of interacting particles, 
a two-body interaction between particles may result in 
the onset of chaos and thermalisation, even if the system 
under consideration consists of a relatively small num- 
ber of particles (see, for example, the reviews || ||, ^] 
and references therein). In application to quantum com- 
puters, quantum chaos may play a destructive role since 
it increases the system sensitivity to external perturba- 
tions. Simple estimates obtained for systems of L inter- 
acting spins show that with an increase of L the chaos 
border decreases, and even a small interaction between 
spins may result in chaotic properties of eigenstates and 
spectrum statistics. On this ground, it was claimed || 
that quantum chaos for a large number of qubits can 
not be avoided, and the idea of a quantum computation 



meets serious problems. 

On the other hand, recent studies || of a realistic 1/2- 
spin model of a quantum computer show that, in the pres- 
ence of a magnetic field gradient, the chaos border is inde- 
pendent on L, and that quantum chaos arises in extreme 
situations only, which are not interesting from the practi- 
cal viewpoint. One should stress that a non-zero gradient 
magnetic field is necessary in the model || for a selec- 
tive excitation of different qubits under time-dependent 
electromagnetic pulses providing a specific quantum pro- 
tocol. 

Another point that should be mentioned in the con- 
text of quantum chaos is that typical statements about 
chaos refer to stationary eigenstates and spectrum statis- 
tics. However, quantum computation is essentially a 
time-dependent problem. Moreover, the time of com- 
putation is restricted by the length of a quantum proto- 
col. Therefore, even if stationary Hamiltonians for single 
pulses reveal chaotic properties, it is still not clear to 
what extent stationary chaos influences the evolution of 
a system subjected to a finite number of pulses. 

In contrast with our previous studies H, in this paper 
we investigate a time evolution of a 1/2-spin quantum 
computer system subjected to a series of pulses. Specifi- 
cally, we consider a quantum protocol that allows to cre- 
ate an entangled state for remote qubits. For this, we 
explore the model in the so-called selective regime, using 
both analytical and numerical approaches. Our analyti- 
cal treatment shows that in this regime there is no any 
fingerprint of quantum chaos. Moreover, we show that 
a kind of perturbation approach provides a complete de- 
scription of the evolution of our system. 

We concentrate our efforts on the introduced quantity 
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(dynamical fidelity). This quantity characterizes the per- 
formance of quantum computation associated with the 
dynamical errors. Dynamical fidelity differs from the fi- 
delity which is associated with the external errors, and 
which is widely used nowadays in different applications 
to quantum computation and quantum chaos, see for in- 
stance JjJ. Our study demonstrates an excellent agree- 
ment of analytical predictions with numerical data. 

The structure of the paper is as follows. In Section II 
we discuss our model and specify the region of parame- 
ters for which our study is performed. In Section III we 
explore a possibility of quantum chaos in the selective 
regime, and analytically show that chaos can not occur 
in this case. We provide all details concerning the quan- 
tum protocol in Section IV, and demonstrate how per- 
turbation theory can be applied to obtain an adequate 
description of the fidelity in dependence on the system 
parameters. Here, we also present numerical data and 
compare them with the predictions based on the per- 
turbative approach. The last Section V summarizes our 
results. 



II. SPIN MODEL OF A QUANTUM 
COMPUTER 



Our model is represented by a one-dimensional chain 
of L identical 1/2-spins placed in an external magnetic 
field. It was first proposed in |8) (see also |)[ If ) 



as a simple model for solid-state quantum computation. 
Some physical constraints are necessary in order to let 
it operate in a quantum computer regime. To provide a 
selective resonant excitation of spins, we assume that the 
time independent part B z = B z (x) of a magnetic field is 
non-uniform along the spin chain. The non-zero gradient 
of the magnetic field provides different Larmor frequen- 
cies for different spins. The angle 9 between the direc- 
tion of the chain and the z-axis satisfies the condition, 
cos 9 = 1/ \fi. In this case the dipole-dipole interaction 
is suppressed, and the main interaction between nuclear 
spins is due to the Ising-like interaction mediated by the 
chemical bonds, as in a liquid state NMR quantum com- 
putation Q. 

In order to realize quantum gates and implement op- 
erations, it is necessary to apply selective pulses to single 
spins. The latter can be distinguished, for instance, by 
imposing constant gradient magnetic field that results in 
the Larmor frequencies u) k = j n B z (xk) = ujQ + ak, where 
7„ is the spin gyromagnetic ratio and ak = Xk and Xk is 
the position of the fc-th spin. If the distance between the 
neighboring nuclear spins is Ax = 0.2 nm, and the fre- 
quency difference between them is A/ = a/2ir = 1 kHz, 
then the corresponding gradient of the magnetic field can 
be estimated as follows, \dB z /dx\ = A//(7„/27r)Ax « 
1.2 x I0 4 T/m. Here we used the gyromagnetic ratio for 
a proton, 7„/27r ss 4.3 x 10 7 Hz/T. Such a magnetic field 
gradient is experimentally achievable, see, for example, 

E HI. 



In our model the spin chain is also subjected to a 
transversal circular polarized magnetic field. Thus, the 
expression for the total magnetic field has the form 

B(t) = [b p L cos(v p t + ip p ),-b p L sm(iy p t + tp p ),B z (x)}. (f) 

As is mentioned above, here B z (x) is the constant mag- 
netic field oriented in the positive x-direction, with a pos- 
itive z-gradient (therefore, a > in the expression for 
the Larmor frequencies). In the above expression, 6^, 
v pi and ip p are the amplitudes, frequencies and phases 
of a circular polarized magnetic field, respectively. The 
latter is given by the sum of p — 1, ...,P rectangular 
time-dependent pulses of length t p+ i — t p , rotating in the 
(x, y)— plane and providing a quantum computer proto- 
col. 

Thus, the quantum Hamiltonian of our system has the 
form, 



L-1 

H = - E 

fc=0 



u k P k + 2 £ J M J£J* 



P L-1 
\ E Qp{t)tt p E e-^*-^/- + e wpt+ifpl+ 
p=l k=0 V 



(2) 



where the "pulse function" <d p (t) equals 1 only during the 
p-th pulse, for t p < t < tp+i, otherwise, it is zero. The 
quantities Jk, n stand for the Ising interaction between 
two qubits , ujk are the frequencies of spin precession in 
the B z — magnetic field, VL p is the Rabi frequency of the 
p-th pulse, I*' V:Z — (l /2)a^ v ' z with tr^' y ' z as the Pauli 
matrices, and if = 1% ± il^. 

For a specific p-th pulse, it is convenient to represent 
the Hamiltonian (||) in the coordinate system that rotates 
with the frequency v p . Therefore, for the time t p < t < 
t p+ i of the p-th pulse our model can be reduced to the 
stationary Hamiltonian, 



k=0 



+ alt - pit) 



2 E Jk,nItHn, 



(3) 

where = (uok — fp), ot = Q p cos(p p , and (3 = Q p simp p . 

We start our considerations with the simplified case 
of the Hamiltonian (3) for a single pulse, by choosing 
ip p = 0. We also assume a constant interaction between 
nearest qubits only ( Jk >n — J5k,k+i), and we put fl p = Q. 
Then the Hamiltonian (3) takes the form, 



L-1 



H <p) = j2 [-n k F k -2Jim + i 



k=0 



L-1 
k=0 



-V. (4) 



In the z-representation the Hamiltonian matrix of size 
2 L is diagonal for Q = 0. For O ^ 0, non-zero off-diagonal 
matrix elements are simply H kn = H n k = —0/2 with 
ji =/= k. The matrix is very sparse, and it has a specific 
structure in the basis reordered according to an increase 
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of the number s. The latter is written in the binary repre- 
sentation, s = ii_2, — , «o (with i s = or 1, depend- 
ing on whether the single-particle state of the i— th qubit 
is the ground state or the excited one). The parameter 
n thus is responsible for a non-diagonal coupling, and we 
hereafter define it as a "perturbation" . 

In our previous studies Q we have analyzed the so- 
called non-selective regime which is defined by the con- 
ditions, Q p Sujk J. This inequality provides the 
simplest way to prepare a homogeneous superposition of 
2 L states needed for the implementation of both Shor 
and Grover algorithms. Our analytical and numerical 
treatment of the model (2) in this regime has shown 
that a constant gradient magnetic field (with non-zero 
value of a) strongly reduces the effects of quantum chaos. 
Namely, the chaos border turns out to be independent on 
the number L of qubits. As a result, for non-selective ex- 
citation quantum chaos can be practically neglected (see 
details in ||). 

Below we consider another important regime called se- 
lective excitation. In this regime each pulse acts selec- 
tively on a chosen qubit, resulting in a resonant transi- 
tion. During the quantum protocol, many such resonant 
transitions take place for different p pulses, with different 
values of v v = ojk ■ The region of parameters for the se- 
lective excitation is specified by the following conditions 



Q P <C Jfc, n <C a <C u>k- 



(5) 



The meaning of these conditions will be discussed in next 
Sections. 



III. ABSENCE OF QUANTUM CHAOS IN THE 
SELECTIVE EXCITATION REGIME 

Here, we consider the properties of the stationary 
Hamiltonian (^) in the regime of selective excitation. In 
order to estimate the critical value of the interaction J, 
above which one can expect chaotic properties of eigen- 
states, one needs to compare the typical value of the off- 
diagonal matrix elements (which is 0/2) with the mean 
energy spacing Sf for unperturbed many-body states that 
are directly coupled by these matrix elements. Therefore, 
the condition for the onset of chaos has the form, 



n 



> s f 



(AE) 



f 



(6) 



Here (AE) f is the maximal difference between the ener- 



gies Eq and Eq 1 ' corresponding to a specific many-body 
state |1) , and all other states |2) of Hq, that have non- 
zero couplings (1|V|2). Correspondingly, Mf is the num- 
ber of many-body states |2) coupled by V to the state 
|1). A further average over all states |1) should be then 
performed. 

In fact, such a comparison (JTJ) is just the perturba- 
tion theory in the case of two-body interaction. Strictly 



speaking, the above condition in a strong sense (fl 6t) 
means that exact eigenstates consist of many unper- 
turbed (V — 0) states. Typically, the components of such 
compound states can be treated as uncorrelated entries, 
thus resulting in a chaotic structure of excited many- 
body states. However, one should note that in specific 
cases when the total Hamiltonian is integrable (or quasi- 
integrable), the components of excited states have strong 
correlations and can not be considered as chaotic, al- 
though the number of components with large amplitudes 
can be extremely large (see details in Q). 

It is relatively easy to estimate Mf in the regime of 
selective excitation. Let us consider an eigenstate of Hq, 
|1, 0,0,0, 1,0, ...,0,0, 1,0), as a collection of 0's and l's 
that correspond to —1/2 and 1/2-spin values. Since the 
perturbation V is a sum of L terms, each of them flipping 
one single spin, one gets Mf ~ L. 

In order to estimate (AE)f, let us first consider the 
action of V on the fc-th spin, and for each spin compute 
the relative energy difference between the final and the 
initial energy. One can find that if the fc-th spin has 
the value 1/2, there are four possible configurations of 
neighbor spins coupled by the perturbation, 



.0,1,0. 
.1,1,1. 
.1,1,0. 
.0,1,1. 



.0,0,0. 
.1,0,1. 
.1,0,0. 
.0,0,1. 



(7) 



If the fc-th spin has the value —1/2, there are also four 
possible different arrangements, 



-o.o.o... -> |... o,i,o...), 

-i.o.i... -> |-i,i,i-), 

-i.o.o- -> |...i,i,o...), 

-(i.o.i- -> |...o,i,i...), 



(8) 



which are the inverse transitions of (0) . Correspondingly, 
the energy changes are determined by the relation, 



\EW-E»>\ = \t k ±2J\,\t k \, k = l,...,L-2. (9) 

The analysis for the border spins can be performed in 
a similar way, and one gets four possible configurations, 
with the following energy changes, 



\E 



(/) 



Ejf>\ = \£ k ±J\, k = 0,L-l. 



(10) 



Summarizing the above findings, and setting for instance 
v p = u)q, one can conclude that (AE)f can be estimated 
as follows, 

(AE) f = Max(\E { f) - E^\) = u L -i - ujq + J- (11) 

As a result, the condition for the onset of quantum 
chaos can be written in the form, 

Q (AE) f _ u £ _! -ojq + J _ a(L-l) + J 

2 > M f ~ L ~ L ' 1 ' 
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or 

f2>fi cr ~a+y. (13) 

However, this critical value is outside the range of pa- 
rameters required to be in the selective excitation regime 
il < a (see inequality (6)). Thus, we can conclude that 
quantum chaos for stationary states can not appear in 
the selective excitation regime. Note that the analysis is 
done for a single pulse of a time-dependent perturbation. 

IV. FIDELITY OF A QUANTUM PROTOCOL 

Analytical results obtained above, show that during a 
single electromagnetic pulse the system can be described 
by perturbation theory. Indeed, if matrix elements of 
perturbation are smaller than the energy spacing between 
directly coupled many-body states, exact eigenstates can 
be obtained by perturbation theory. Thus, one can ex- 
pect that for a series of time-depended pulses the evolu- 
tion of the systems can be also treated making use of a 
kind of perturbative approach. 

In what follows, we study the dynamics of the system 
by applying a specific set of pulses (quantum protocol) 
in order to create an entangled state for remote qubits 
(with k = and k = L — 1) starting from the ground 
state, |V>o) = |0l-i, Oi, Oo) (we omit the subscripts 
below). Our main interest is in estimating the errors that 
appear due to unwanted excitations of qubits. We show 
that these errors can be well understood and estimated 
on the basis of the perturbation theory developed for our 
time-dependent Hamiltonian (j^) , in the parameter range 
where the protocol holds. 

A. Selective excitation regime and perturbation 
theory 

Any protocol is a sequence of unitary transformations 
applied to some initial state in order to obtain a final 
state, IV' 1 ). In this model of quantum computer the 
protocol is realized by applying a number of specific rf- 
pulses, so that we get a state \ip r ) which is, in principle, 
different from the ideal state \ip 1 )- The difference be- 
tween the final state \i[i r ) and the ideal state can be 
characterized by a dynamical fidelity, 

F=\mr)\ 2 . (w) 

Note that, in our case the dynamical fidelity F does not 
explicitly depend on a perturbation parameter added in 
the Hamiltonian (2) in order to get a distorted evolution, 
as is typically assumed in the study of quantum chaos. 
Indeed, the final state is determined by the total Hamil- 
tonian (2), 

JL A _j f tp H(t)dt 

\r) = U(T)\^ ) = l[Te Jt »- |^ ), (15) 
P =i 



where T = t p is the total time to entangle spins, U(T) 
is the unitary operator given by the sequence of pulses 
in the protocol, and T— is the usual time-ordered prod- 
uct. Therefore, it is not possible to identify a single per- 
turbation parameter which is responsible for a "wrong" 
evolution of the system. 

The selective excitation regime is characterized by the 
action of pulses that are resonant with a transition be- 
tween two energy states which differ for the state (up 
or down) of one spin only. A close inspection of the 
time independent Hamiltonian (^), defines the region of 
parameters where the selective excitation of single spins 
can be performed. 

Diagonal elements of Hamiltonian (0) are given by the 

eigenvalues Eq of Hq, while non-zero off-diagonal ele- 
ments are constant and equal to —0/2. In order to have 
a resonant transition between two energy states, their en- 
ergy difference A has to be zero. On the other side, for 
each state no more than one resonant transition should 
be allowed. So, we require that the energy differences 
given by Eqs. (p|,^0|) are different from zero (apart from 
the wanted transition) . This leads to the following set of 
equalities ("fake transitions"), 



J 


=1 

= 4 


when 


k = 


1,. 


..,L- 


3, 


J 


when 


k = 


1,. 


..,L- 


3, 


J 


= ak 


when 


k = 


1,.. 


,.,L- 


2. 


J 


= 4 


when 


k = 


1,. 


..,L- 


2. 



From Eqs.(fl6|) it is easy to see that first "fake" transi- 
tion appears for Ji = a/4, the second for J 2 = a/2 and 
so on up to the last one for Jf = a(L — 2). All these 
resonances can be avoided if we choose a ^> 4 J. Indeed, 
it is not enough to chose simply a > 4 J since resonances 
have finite widths. 

Transitions can be defined according to their energy 
difference A, 

1) resonant transitions, A = 0; 

2) near-resonant transitions, A ~ J; 

3) non-resonant transitions, A ~ a. 

For a 3> 4J, each state can undergo one resonant or 
near-resonant transition only, and many non-resonant 
ones. The latter can be neglected if we choose a 3> 
il. Under these conditions we can form couples of 
states, connected by resonant or near-resonant transi- 
tions, and we can rearrange the Hamiltonian matrix (Q) 
by 2 x 2 block matrices representing all resonant and 
near-resonant transitions. This allow us to describe the 
dynamical evolution of the system as a two-state prob- 
lem. 

Using this procedure, the entire sequence of pulses can 
be evaluated. Note that special attention has to be paid 
to an additional phase shift that arises between any two 
pulses, due to the change of frame. We remind that the 
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transformation between the rotating and the laboratory 
frame is given by the expression, 



m))Lab = 



' k Mt)) 



Rot 



(17) 



Indeed, let us consider an initial basis state \m) at 
time t = 0, and find the probability for a resonant (A — 
0) or near-resonant (A ~ J) transition to the state \p) 
with the energy difference E p — E m . Here E p and E m 
are the eigenenergies of the time-independent part of the 
Hamiltonian (^), written in the laboratory frame. 

Setting, 



m = £ Cn(*)e-* Bnt |n), 



(18) 



and c p (0) = 0, the application of a pulse for a time r 
gets: 



n(r) = c m (0)[cos(4f) + i f 
>(r) = c m (0)[i § sin(^)]e i 



sin(f )]e" 



-iE„ 



^ (19) 

where A = \/VL 2 + A 2 . Note that Eqs.(jl|) refer to the 
laboratory frame. 

As we can see, the parameter e determined as 



(20) 



characterizes the probability of resonant and near- 
resonant transitions. In particular, the probability of un- 
wanted near-resonant transitions goes like e, and it can 
be reduced by assuming J 3> 51. Combining all the above 
expression, we get the condition (|J). 

Correspondingly, the probability for a non-resonant 
transition (neglecting terms of the order 1 / L, and as- 
suming a 3> SI ) is given by the parameter r\ |ll[ , 



n 



4a 2 ' 



(21) 



We would like to stress that even if the ideal state has 
been constructed taking into account resonant transitions 
only, our dynamical fidelity is a measure of dynamical 
errors that are due to near and non-resonant transitions. 

Let us now briefly discuss the perturbative approach 
that is based on recent studies published in Ref. |ll|] . The 
main idea is that for each p-th. pulse the unperturbed ba- 
sis can be rearranged in such a way that the Hamiltonian 
matrix is represented by 2 x 2 block matrices, as described 
above. This is what we call unperturbed Hamiltonian for a 
specific p-th pulse. One should note that this unperturbed 
Hamiltonian is Sl-dependcnt. Let us now define by V the 
Jl-dependent part which is responsible for non-resonant 
transition and not described by the 2x2 block matrices. 
Then it is easy to obtain the unperturbed eigenstates, 
\ipq), and the unperturbed eigenvalues, e°, by diagonal- 
izing each of the 2x2 blocks independently. 



After this step, one can compute the perturbed eigen- 
states by taking into account first order terms only, 



,n ,o Wa'h 



(22) 



Note that this perturbative approach is supposed to 
be valid when Eqs.([l6]) are not satisfied, and when the 
errors due to near-resonant transitions are much larger 
than the errors due to non-resonant ones, e 3> r\. 



B. Quantum protocol 

Let us briefly sketch the algorithm and the particular 
protocol which was developed in Ref. [jl(| . Starting from 
the ground state |V>o) = 1 0. . .0) and applying a number 
of specific pulses, we would like to generate the following 
entangled state, 



1 



|y/) = — (|o...o) + |io. 



.01)). 



(23) 



This algorithm could serve, for instance, as the first step 
for a more general teleportation protocol, and for an im- 
plementation of conditional quantum logic operations. 

The algorithm can be realized in the following way (for 
details see Ref. [fl0|), 



|0,...,0) -(|0,...,0> + |1,0,..,0» 
-(|0,...,0> + |1,1,0,..,0» 
^(|0,...,0) + |1,1,1,0,..,0)) 
^(|0,...,0) + |1,0,1,0,..,0)) 
-»...-» (|0,...,0) + |1,0,.. I 1)). 



(24) 



Physically, the above algorithm can be done by apply- 
ing suitable r/-pulses that are resonant to the desired 
transitions. The latter are originated from induced Rabi 
oscillations between the resonant states. 

To flip the fc-th spin we have to choose the frequency 
v of the rf-pulse according to the relation v p — E\ — E% , 
where |1) and \2) are the states involved in the transi- 
tion and E\ and Ei are the eigenenergies of the time- 
independent part of Hamiltonian (||). For instance, for 
the first pulse we have to set, v\ — \E\i q q) ™ ^|o,...,o)l> 
and we have to apply it for a time t\ = tt/20, to get 
equal superposition of the states involved in the tran- 
sition. For other pulses we require that the first state 
(|0,...,0)) remains the same (apart from an additional 
phase), while the second state flips the k— th spin. This 
amounts to say that the probability of unwanted states 
is due to non-resonant transitions of both states, and to 
near-resonant ones only of the first state in the r.h.s. of 
Eq. (p4[) . Specifically, the state |0,...,0) undergoes near- 
resonant transitions with A = 2 J for each pulse, except 
the first one which is resonant, and the forth for which 
A = 4 J. One should also take into account that at each 
pulse the state |0, 0) get an additional phase, see eqs. 
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( |l9| ) . We took them into account in the definition of the 
ideal state. 

Since in the selective excitation regime we have e 3> 
v, contributions from near-resonant transitions are much 
larger than the one due to non-resonant transitions. Our 
algorithm consists of 2L — 2 separate pulses, therefore, 
some modifications are necessary in order to be able to 
control small unwanted probability. For the product of 
probabilities this requires to impose, 2Le < 1 and 2Lr\ < 
1, that can be written in the form, 



n 
7 



< 



2 

V 



< 



(25) 



where L ^> 1. 

Before discussing our numerical results we would like 
to stress that in contrast to what is mainly considered 
in the literature, the time for our dynamical fidelity is 
not an independent variable. Indeed, the length of the 
protocol is determined by the total number of qubits, 
L. Specifically, 2L — 2 pulses are necessary in order to 
create the entangled state, so that the protocol time T is 
proportional to the number of qubits. 



C. Dynamical fidelity: theory and numerical data 

Quite unexpectedly, the dynamical fidelity (14) in- 
creases with an increase of the Ising coupling J, as soon 
as J < a/4. This is due to the fact that the probability 
of unwanted near-resonant transitions is proportional to 
e ~ (fl/J) 2 (see Eq.(pO|) where A ~ J for near-resonant 
transitions ) . The larger is J, the smaller is the probabil- 
ity of non-resonant transitions, therefore, the dynamical 
fidelity F increases with an increase of J. 

In Fig.[j] we show how the dynamical fidelity (|14| ) de- 
pends on the inter-qubit interaction J. For convenience, 
the function 1 — F is shown here and below, instead of 
F. Numerical data have been obtained in two different 
ways. Full curve corresponds to exact computation of 
the time-dependent Hamiltonian (||). Data in Fig.^Ja are 
compared with those obtained from the perturbative ap- 
proach explained above. 

Apart from very strong peaks (see Fig.|l|b) for which 
the dynamical fidelity vanishes, one can say that the 
global tendency is an improvement of the dynamical fi- 
delity for larger values of J. However, strong oscillations 
occur reflecting a resonant nature of the dynamics of our 
system. Perfect agreement between perturbative results 
and numerical data is found for very large variations of 
the interaction strength J. 

High peaks for 1 — F, clearly seen in Fig.|l|b, are due 
to Eqs.([l6|), in these points the quantum algorithm fails. 
Thus, one should avoid these situations in a quantum 
computation. As for the minima in Fig.|l|a for which the 
dynamical fidelity is close to one, they occur when e = 0, 
or, when 
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FIG. 1: The dependence 1 — F is shown as a function of the 
Ising coupling J for L = 6 spins, Q — 0.118 and a — 100. Full 
line represents the numerical data for the dynamical fidelity 
F defined by Eq.(jl^), and obtained from direct numerical 
computation of the system evolution, (a) Full circles stand 
for perturbative calculations, and full curve corresponds to 
numerical results, (b) The same numerical results as in (a), 
but for a larger range of J. The theoretical expression as 
given by Eq.(|30j) is also shown in (a). 



where k is an integer number. This relation corresponds 
to the 27rfc-condition ju| H), 0]. 

Let us now explore the dependence of the dynami- 
cal fidelity on the parameter a which is proportional 
to the gradient of the external magnetic field, a = 
'y n Ax(dB z (x)/dx), where Ax is the distance between 
neighboring qubits. (Below, we shall refer to the param- 
eter a as to the magnetic field gradient.) Numerical data 
for the dependence of 1 — F on a are presented in Figj^. 
One can see that the dynamical fidelity is getting bet- 
ter for large enough values of a. We already mentioned 
that for a < 4 J a problem may arise in the protocol due 
to "fake" transitions. On the other side, in the regime 
a 4J the dynamical fidelity reaches an asymptotic 
value which depends on J and £1 only, see Fig.y. 

It is also important to understand the dependence of 
the dynamical fidelity on the Rabi frequency f2. The data 
manifest two specific properties demonstrated in Figj^. 
The first one is a global decrease of the dynamical fidelity 
with an increase of f2. The second peculiarity is due to 
strong oscillations that occur for e = 0, namely for those 
Q values that correspond to the 27rfc-conditions, 
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For these values of near-resonant transitions vanish, 
and non-resonant transitions remain only. Thus, the dy- 
namical fidelity has maxima which provides, in principle, 
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FIG. 2: The dependence 1 — F as a function of the gradient 
of a magnetic field is shown for L — 6 spins, Q = 0.118 and 
J — 1. As one can see, for a < 4J the correspondence is not 
so good as for a > 4J. 
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FIG. 3: The difference 1 — F as a function of the Rabi fre- 
quency Q for L — 6 spins, J = 1 and a = 100. Full curve is 
the result of direct numerical simulation, circles are obtained 
from the perturbative approach described in the text. Arrows 
show few resonant values of f2 given by Eq.(|26|). 



the best condition for a quantum computation. 

Nevertheless, let us consider the values of fl that cor- 
respond to maxima in Fig.(^). This we do in order to 
make an estimate in the worst possible condition. A 
brief analysis of the fidelity for specific values fl = ttf. 
will be sketched in the last subsection. As one can see, 
for values of different from Qfe, the "average" dynami- 
cal fidelity increases when the Rabi frequency decreases. 
This is due to the fact that the probability to gener- 
ate unwanted states (due to both non-resonant and near- 
resonant transitions), is proportional to (il/A) 2 . There- 
fore, the smaller is the more reliable the algorithm is. 
Note that the agreement with the perturbative approach 



is excellent. 

On the other hand, we cannot choose an extremely 
small value of fl since it implies a large time duration of 
the pulse (t ~ ir/£l). Note, that the total time for a quan- 
tum protocol should be kept well below the decoherence 
time (the latter can be quite large for nuclear spins |jT5f ). 
Taking this point into account, an optimal choice is to 
take the largest possible value, f2 = O2 = 2 J/ \/T5 < J, 
and large enough value of a (in order to significantly sup- 
press the non- resonant transitions). 



D. Fidelity: dependence on the number of qubits 

Finally, we studied the dependence of the dynamical 
fidelity on the number L of spins in the chain. As was 
noted before, for a chosen protocol its length is propor- 
tional to L. Numerical data clearly manifest a linear 
decrease of the dynamical fidelity with the number of 
qubits, see Fig. |[ 

This dependence can be also understood on the basis of 
perturbation theory, by neglecting the change of phases 
between the ideal and the real state. This can be done 
since we already considered the main change of phases in 
the definition of the ideal state. Let us write, 
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and assume that after M pulses the components Cq and 
ci are, 
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Here cq and ci are the complex amplitudes of two states 
in Eq. (|2^) , but obtained as the result of the quantum pro- 
tocol. This means that the probability of near-resonant 
states after M pulses is small, Me <C 1. Then we get, 
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(2 - Me + 2Vl - Me) ~ 1 - M- 



(29) 



where e ~ f2 2 /4 J 2 and M = 2L — 3 for our specific algo- 
rithm. Combining all the above expressions, one gets, 
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(30) 



which implies a linear decrease of the dynamical fidelity 
with an increase of the number of qubits, L. The slope is 
given by the parameter m t h = — fl 2 /4:J 2 . Of course, this 
consideration is valid far from the resonant conditions 
©• 

To compare numerical data with the above analytical 
estimate for mth, we obtained the slopes in Fig. f|by the 
best fit to a linear dependence. These slopes versus J are 
shown in Fig. |5[ together with the theoretical result for 
m t h- 

As one can see, the agreement is very good except for 
small values of the Ising interaction, J ~ f2, where the 
probability of near-resonant transitions becomes large. 
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FIG. 4: The dynamical fidelity as a function of the number 
L of spins, for different J values and Q — 0.118, a — 100. 
Numerical data (triangles for J = 1.945, circles for J = 5.01 
and squares for J — 9.99) are compared with the results from 
the perturbation theory (crosses). Also shown are the best 
linear fits (dot-dashed lines). 
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FIG. 5: Comparison between theoretical and numerical lin- 
ear slopes for the fidelity, as a function of the interaction J. 



E. Optimal algorithm 

Choosing fl values as given by Eq. (prj|) , one gets that 
the probability for near-resonant transition is zero, e = 
0. So, only non-resonant transitions lead to unwanted 
states. In Fig.(^) we show the fidelity as a function of 
the number of spins L for ilk — 0.1216. These data 
should be compared with the analogous ones indicated 
by triangles in Fig.([I|). 

As one can see, despite the closeness of these two fl 
values (less than 3% of difference), the fidelity is better 
by two order of magnitude (see different scales on the y- 
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FIG. 6: The dynamical fidelity as a function of the number 
L of spins, for J = 1.945 and ft fc = 0.1216, a = 100. Also 
shown is the best linear fit with the slope 9.2 ± 0.3 x 10 -6 . 



axis). It is clear that such preferred O values should be 
chosen in any practical implementation of the algorithm. 
However, due to a high instability of such resonant val- 
ues, see Fig.(||), the detailed analysis can only be done 
within a more general study under the presence of small 
variations in parameters such as f2, J, a. This study is 
currently in progress. 



V. CONCLUSIONS 

We have studied the model of a quantum computer 
that consists of a one-dimensional chain of 1/2-spins 
(qubits), placed in a time-dependent electromagnetic 
field. The latter is given by a sequence of rf-pulses, cor- 
responding to a chosen quantum protocol that allows to 
generate an entangled state for remote qubits from the 
initial ground state. Main attention is paid to the analy- 
sis of the dynamical fidelity, defined as an overlap of the 
actual finite state, with the ideal one determined by the 
quantum protocol. 

We considered the region of the selective excitation 
where the resonant excitations of specific qubits can be 
implemented by time-dependent pulses. Analytical treat- 
ment of the stationary Hamiltonian which describes the 
evolution of the system during a single pulse, has revealed 
that in the selective regime the quantum chaos can not 
appear. Moreover, in this regime a perturbation theory 
can be applied to all quantities of interest. 

Our detailed study of the dynamical fidelity manifests 
excellent agreement between numerical data and the pre- 
dictions obtained in the perturbative approach. In par- 
ticular, we have found how to choose parameters of the 
model in order to get the best dynamical fidelity for the 
creation of the remote entangled state. Specific attention 
has been paid to the dependence of the dynamical fidelity 
on the number L of qubits. We show, both analytically 



9 



and numerically, that the dynamical fidelity decreases 
linearly with an increase of L, and we give an analytical 
estimate for the slope of this dependence. 
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